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We systematically compare filtering methods used to extract topological structures on 577(3) lat- 
tice configurations. We show that there is a strong correlation of the topological charge densities 
obtained by APE and Stout smearing. To get rid of artifacts of these methods, we analyze struc- 
tures that are also seen by Laplace filtering. This combined analysis shows that the topological 
charge density is more fragmented in the presence of dynamical quarks. 
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1. Filtering Methods and Topological Charge Density 

Many methods have been developed to extract the IR content from lattice data. Unfortunately, 
all these methods introduce ambiguities and parameters. Thus, to get a coherent picture of the topo- 
logical structure of the QCD vacuum, it is necessary to find ways of controlling or even removing 
these ambiguities. 

One of the first attempts to filter out the UV 'noise' has been APE smearing [1], defined as: 
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where a determines the weight of the old link and the sum of the attached staples. Throughout this 
work we have used the standard value a = 0.45. The right hand side has to be projected back to the 
gauge group. The projection onto SU (3) is not unique, but a common choice is to define Psu(3) (W) 
as the that element of the group V G St/ (3) that maximizes ReTrjVW^} . Stout smearing [2] 
avoids this projection by using the exponential map U^ tout = exp{z Q^(U,p)} ■ U^ ld , where 
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is a hermitian traceless matrix constructed from all plaquettes containing the old link U^ d and a 
smearing parameter p and hence e'@ G SU (3). 

Both smearing methods are iterative procedures, which lead to smoother configurations in each 
step. But it is hard to decide how many smearing steps one can apply without loosing information 
about the underlying structure. 

A relatively new method is Laplace filtering [3]. This method is based on a truncated spectral 
decomposition of the links in terms of eigenmodes <J>„ (x) of the covariant lattice Laplacian: 
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As the eigenvalues are related to some energy scale (squared), this procedure acts as a low-pass 
filter in the sense of a Fourier decomposition: with lower values ,/V of included modes the Laplace 
filtering gets stronger. At this point it should be stressed that Laplace filtering is completely differ- 
ent from smearing, because it is based on extended objects, namely the eigenmodes, and does not 
locally modify the gauge links in contrast to smearing. 

In this sense this method is quite similar to Dirac filtering [4] . However, Laplacian modes are 
chirality blind and do not obey an index theorem, which connects the zero-modes of the (chiral) 
Dirac operator to the topological charge. From a computational point of view Laplace filtering is 
much cheaper than Dirac filtering (with good chirality properties). 



'The maximum is found iteratively using code included in the CHROMA software package for lattice QCD [5]. 
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lat. size 


lat. spacing 




quenched 


16 3 -32 


0.148 


7.90 


dynamical 


16 3 -32 


0.150 


4.65 



Table 1: Details on the gauge configurations used in this work. 

After filtering we measure the topological charge density q(x) = (l/\6n 2 )Tr{F llv (x)F ,lv (x)}, 
where F„ v is an improved discretization of the field strength tensor [6], which combines 1x1,2x2 
and 3x3 loops to achieve 0(a) -improvement at tree-level 2 . 

In order to get results for dynamical and quenched configurations which can be easily com- 
pared, we have chosen two SU(3) ensembles with the same lattice spacing and the same physical 
volume (see Tab. 1). The ensembles were generated with the Liischer-Weisz gauge action and a 
chirally improved Dirac operator. For the dynamical simulations two flavors of mass degenerate 
light quarks were used (details can be found in [7]). 



2. Comparison of the Methods 

For quenched SU (2) lattice configurations we had observed a similarity of the topological 
density obtained with the different filtering methods [8]. For a quantitative comparison one needs 
a measure of the similarity. In [8] the following quantity was introduced: 

Z AB ^^^- (2.1) 
Xaa Xbb 

with Xab = (1/^)L,y(<7a(X) — ^aX^bM ~~ <7s) being the correlator of two topological charge den- 
sities q A (x) and qs{x) (mean values q = Q/V are subtracted for convenience). 

E is a positive quantity and equals 1 , if the densities differ only by a constant scaling factor. 
The great advantage of this quantity is that we do not have to know the normalization factors of the 
topological charge densities, since these factors drop out in eq. (2.1) and we are able to compare 
only the relative difference at each lattice point. 

Moreover, we can define 'best matching' pairs of filter parameters, by maximizing E A b- To 
be precise, we keep the smearing parameters constant and compute E A b for up to 50 APEsteps and 
50 Stout steps respectively up to 500 Laplace modes. More smearing steps have not been applied, 
as too much smearing will destroy the local topological structure, whereas more Laplace modes 
would be too expensive. 

The matching of APE and Stout smearing, see Fig. 1(a), shows a one-to-one correspondence 
of the topological charge densities [9, 10]. Thus, the different smearing techniques reveal almost 
the same structures, but they will also find the same unphysical artifacts. To control the ambiguities 
it is necessary to complement smearing by an independent filtering method, in our case Laplace 
filtering. 

The comparison is not so perfect when we consider APE smearing and Laplace filtering, see 
Fig. 1(b). Also here we find some kind of ridge line of the best matching values of S. This ridge 



for details on the improvement coefficients see [6] 
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(a) APE vs. Stout (b) APE vs. Laplace 



Figure 1: Comparison of APE with Stout smearing and APE smearing with Laplace filtering for an ensemble 
dynamical configurations. The degree of similarity S (see the definition (2.1)) is given by the contour lines. 
Shown are mean values of 80 configurations for APE vs. Stout smearing and 30 configurations for APE 
smearing vs. Laplace filtering. We find for APE vs. Stout smearing a band with E > 0.95. Thus, both 
methods find the same structures and artifacts. The values of S for APE smearing vs. Laplace filtering does 
not reach such high values. For mild filtering we have plotted a ridge line of best matching parameters that 
will be used in the cluster analysis later. The described behavior applies also to quenched configurations . 

is quite pronounced for weak filtering, but we are confronted with the problem to identify the best 
matching parameters in the strong filtering regime, i.e. for less than 50 Laplace modes and more 
than 20 smearing steps. Furthermore, there is no set of best matching parameters with S > 0.5 in 
this regime and thus one can hardly say that both methods agree at all (E = 0.5 roughly means that 
only every second object is found by both filters). This finding flags a warning that any kind of 
smearing should only be applied with great care. 

3. Cluster Analysis of the Topological Charge Density 

In order to obtain information about the topological structure of the QCD vacuum that can be 
compared to continuum models, we analyze the cluster structure of the topological charge density. 
Two lattice points belong to the same cluster, if they are nearest neighbors and have the same sign 
of the topological charge density. 

For such clusters we had found an interesting power law [8, 10]. To that end one cuts the 
absolute value of the topological charge at a (variable) cut-off q cut and considers the number of 
clusters above this cut-off as a function of the fraction / of all points obeying \q(x)\ > q cut w.r.t. 
the total number of lattice points. 

The exponent £, = dlogA r c i ust (^ cut )/dlogA'p i nts (<7 cu t) of this power law is highly characteristic 
for the topological structure of the QCD vacuum. Different models lead to different predictions. 
This allows for a very sensitive test. 

Pure noise, for instance, yields an exponent £ = 1, because every point forms its own cluster. 
On the other hand, the exponent is close to zero for very smooth densities with large structures, 
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Figure 2: Exponent | for the clusters found by APE smearing (a), Laplace filtering (b) and for a matched 
filtering of both methods (c). The exponents predicted by the dilute instanton gas have been included in all 
plots. In the latter plot the measured exponents have a step-like behavior, because the number of smearing 
steps is an integer quantity and we have to match several numbers of Laplace modes to the same number of 
smearing steps. The smooth lines are interpolations corresponding to a "continuous" smearing parameter. 
The number of smearing steps (in parentheses) refers to the matching of quenched configurations and is 
slightly different for the dynamical case. Errors from the ensemble average have been included, but are 
partly too small to be seen. 



because in such a situation a small change in the cut-off will add new points but not reveal new 
clusters. 

The exponent can be related to the size distribution d(p) ~ p@ of equally charged topological 
objects with arbitrary shape function and arbitrary dimensionality [8]. For instance a dilute gas of 
instantons would have B, = 0.64 in the quenched case and % = 0.66 in the dynamical case (Nf = 2). 

A cluster analysis has one great advantage. It allows to reduce ambiguities coming from a 
single filter by taking only those clusters into account, which are common to different filters. If 
there is an artifact coming from one method, it is unlikely that this artifact will also be seen by the 
other one, such that the common structures are almost free of ambiguities. 

The exponent for clusters common to APE and Stout smearing can be found in Fig. 2(a). 
Obviously, the exponents of the dynamical configurations lie above the quenched values. The 
difference of the cluster exponents quenched vs. dynamical vanishes for stronger smearing (~ 30 
steps) and the exponents settle down to almost the same plateau. So we have reasons to believe that 
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Figure 3: Total number of distinct clusters for a constant fraction / = 0.0755 of points lying above the 
cut-off with smooth interpolation. 



too much smearing destroys the impact of dynamical quarks. 

In Fig. 2(b) the same analysis is done for Laplace filtering (only). We find a slightly bigger 
difference of the exponents of dynamical and quenched configurations. Furthermore, this differ- 
ence remains, unlike for APE smearing, at every stage of filtering. At first sight the value for 50 
modes seems to be odd, but this is just an artifact of the incomplete reconstruction of the topologi- 
cal background. For a small number of modes, we get a very spiky structure and we will find many 
new clusters, if we lower the cut-off. Thus, we get a higher exponent. 

To get rid of the ambiguities of APE smearing and Laplace filtering, we perform in Fig. 2(c) a 
matched cluster analysis, i.e. we consider only those clusters, which are common to both methods. 
We use the optimal set of filter parameters according to the maximal values of S (and corresponding 
to points on the ridge line of Fig. 1) on the abscissa. 

Fig. 2(c) shows that the cluster exponent is larger in the dynamical case than in the quenched 
case. According to the considerations from above, this implies a larger exponent /3 of the size 
distribution. Hence, very small topological objects become suppressed when quarks are taken into 
account. 

In the plots we have included the values of the exponent in instanton gases. These have the 
same ordering with a slightly smaller difference. Generally our measured values of the exponent 
are not far off the instanton gas values (which was not the case in SU (2) [8]). 

Another interesting quantity is the absolute number of clusters above a certain cut-off, see 
Fig. 3. In the dynamical case we find almost twice as many clusters. Hence, we observe a more 
fragmented topological structure in the presence of dynamical quarks. 

Our findings can be compared to the results obtained by the Adelaide group [11]. They have 
observed an "increasing density of nontrivial field configurations" and a suppression of small in- 
stantons. Using a completely different approach, we come to the same conclusion. The advantages 
of our analysis are that our dynamical fermions have nicer chiral properties than the staggered ones. 
Moreover, we do not have to postulate any shape function of the topological objects. Therefore, 
our method works for all topological building blocks and with our matched analysis we can obtain 
results which are almost free from ambiguities of the filtering process. 
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4. Conclusions 

We have performed a comparative study on the local topological structure of quenched and 
dynamical configurations seen by various filtering methods, namely APE smearing, Stout smearing 
and Laplace filtering. 

To this end we have analyzed the power-law behavior of topological charge clusters seen by the 
individual methods and common to both methods. This analysis shows clearly a larger exponent in 
the presence of dynamical quarks, which can be related to a larger exponent of the size distribution 
d(p) ~ p' 3 of topological objects. Hence, very small topological objects become suppressed in 
the dynamical case. Moreover, the total number of clusters for a constant total cluster volume is 
substantially higher for dynamical configurations. Thus the topological structure of the vacuum is 
more fragmented in the presence of fermions. 
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